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ABSTRACT 

We present the first study of the formation and dissipation of current sheets at electron scales 
in a wave-driven, weakly coUisional, 3D kinetic turbulence simulation. We investigate the rela- 
tive importance of dissipation associated with collisionless damping via resonant wave-particle 
interactions versus dissipation in small-scale current sheets in weakly collisional plasma turbu- 
lence. Current sheets form self-consistently from the wave-driven turbulence, and their filling 
fraction is well correlated to the electron heating rate. However, the weakly collisional nature 
of the simulation necessarily implies that the current sheets are not significantly dissipated via 
Ohmic dissipation. Rather, collisionless damping via the Landau resonance with the electrons is 
sufficient to account for the measured heating as a function of scale in the simulation, without the 
need for significant Ohmic dissipation. This finding suggests the possibility that the dissipation 
of the current sheets is governed by resonant wave-particle interactions and that the locations of 
current sheets correspond spatially to regions of enhanced heating. 

Subject fieadings: turbulence — plasmas — solar wind 



1. Introduction 

Turbulence plays an important role in space 
and astrophysical plasmas by mediating the trans- 
fer of energy from large-scale motions to the small 
scales at which the turbulence can be dissipated. 
A major unsolved problem is the identification 
of the physical mechanisms that dissipate the 
small-scale turbulent motions, ultimately convert- 
ing the turbulent energy to plasma heat. The 
dynamics at the dissipative scales are typically 
weakly collisional in diffuse astrophysical plas- 
mas, such as the solar wind, so the mechanisms 
responsible for the dissipation and plasma heat- 
ing are described by kinetic plasma physics. Two 
mechanisms have been proposed to be the dom- 
inantly involved in the dissipation process for 
plasma turbulence: collisionless wave-particle in- 



terac t ions (Howes et al.ll2008aHSchekochihin et al 



20091 : iTenBarge et al.l l2013a ) and dissipat i on in 
small-sc ale current sh e ets (IDmitruk et al. _ 2004t 
Markov skii fc Vasauej 12011 : iMatthaeus fc Velli 



20111: lOsman et al.l 120111 



Karimabadi et al.ll2013[l . 



Servidio et al 



20111 



In weakly collisional plasmas, it is well known 
that wave-particle interactions lead to signifi- 
cant collisionless damping of the linear kinetic 
wave modes via the Landau and cyclotron res- 
onances. In turbulent astrophysical plasmas, it 
has been proposed that the fluctuations at per- 
pendicular scales smaller than the ion Larmor 
radius, k±pi > 1, have properties typi ca l of ki 
nctic Alfvcn waves ("Howes et al.' '2008b''a'; 'Howe 
2008; Schckochihin et al. 2009; Salem ct al. 201 
and will therefore suffer collisionless damping 
In the case of the solar wind, it has been sug- 



Rested (IHowes et alj 2008ai ISchekochihin et al . 
2OO9I : iHowes et al.ll2011ai bl) that electron Landau 
damping dominates the dissipation of these tur- 
bulent electromagnetic fluctuations at k^pi > 1. 
Free energy transferred conservatively to the par- 
ticle distribution functions by wave-particle in- 
teractions is ultimately thermalized by arbitrarily 
weak collisions through the action of an entrop y 
cascade in phase space |Schekochihin et al.ll200S) . 
Adopting this paradigm, predictive cascade mod- 
els of the plasma turbulence and its dissipation can 
be constructed by assuming that the coUisionless 
damping rates of the turbulent electromagnetic 
fluctuations are relat ed to the lin ear dam ping rat es 
from kinetic theory (JHowes et a l. 2008a. l2011bl ). 

On the other hand, a number of recent studies 
focusing on the intermittent structures that inher- 
ently develop in plasma turbulence have suggested 
that dissipation dominantly occurs in coherent 
structures, in particular, small-scale current sheets 
([Pmitruk et al. 2004; Markovskii fc Vasquez 2011 



Matthaeus fc Velli 2011: Osman et al._|2011a: 



Servidio et al.ll201ltlKarimabadi et al.ll2013[ ). The 
kinetic physical mechanism by which dissipa- 
tion occurs in current sheets has not been elu- 
cidated. Hybrid kinetic-ion and fluid-electron 
simulations in 2D sugges t perpendicular ion 



heating in current sheets ( Parashar et al.l l2009t 



Markovskii k Vasquedboill ). 2D and 3D Particle 
in-Cell (PIC) simulations of reconnection sug- 
gest the acc eleration of elect r ons by parallel elec- 
tric fi elds dPritc hctt' '2006'; 'Eged al et all l2008l 
2OO9I : iLe et al.l I2OIO: Egedal et all I2OI0I mm . 
and 2D gyrokinetic simulations suggest linear 



phase mixing/Landau damping (iLoureiro et al 
20131 ). Temperature measurements in the near- 



Earth solar wind ha,v e been used to both support 
(lOsman et al.l2011allbl ') and refute (JBorovskv fc Denton 



20111 ) the proposal that plasma heating domi- 
nantly occurs in current sheets. 

In this Letter, we present a wave-driven, 3D gy- 
rokinetic turbulence simulation at scales smaller 
than the ion Larmor radius that self-consistently 
generates small-scale current sheets. We find that 
the current sheet filling fraction is well correlated 
with the electron heating rate in the simulation. 
Yet, Ohmic dissipation is negligible, and the mea- 
sured electron heating rate by scale is well repro- 
duced by assuming dissipation is entirely associ- 
ated with electron Landau damping, suggesting 



the possibility that the current sheets are damped 
collisionlessly by resonant wave-particle interac- 
tions and correspond to regions of local heating. 

2. Kinetic Simulation 

The simulation was performed wi th the Astro- 
physi cal Gyrokinetics code, AstroGK (JNumata et al, 
2OIOI), a triply periodic Eulerian slab code that 
solves the five-dimensional, electromagnetic gy- 
roaveraged Vlasov-Maxwell equations for the 
perturbed gyroaveraged distribution for each 
specie s and the gyroaveraged Maxwell's equa - 
tions (jPrieman fc Chen 1982: Howes et al.ll2006[ ). 
Collisions are treated using a fully conservative, 
linearized, and gyroaveraged collisi on operator 
(JAbel et al.ll2008HBarnes et al.ll2009l) . The simu- 
lations are driven at the simulation domain scale 
with an oscillating Langevin antenna coupled to 
the component of the vector potential parallel 
to the equilibrium ma gnetic field, Bq = -Bqz 
(JTenBarge et al.ll2013b[ ). 

The simulation models a proton and electron 
plasma with a realistic mass ratio mi/m^ ~ 1836, 
/?i = 1, and Ti/Tf, = 1, where Pi = Vu/ v\, va 
is the Alfven speed, and vu — \j2Tijrai is the 
ion thermal speed. The simulation employs a pe- 
riodic domain of size L\ x i^, elongated along 
the straight, uniform equilibrium magnetic field 
Bq. Relevant parameters for the simulation are 
k^Pi e [5,105], k^pije e [1,16], A^ftpiB^ = 
0.2, v^pi/vue = 0.2, and v^pijvue = 0.5, where 
Pi = ^ti/^i is the ion Larmor radius, J7,j is the 
ion gyrofrequency, e — 2TTpi/Lz — 0.081 <; 1 is 
the gyrokinetic expansion parameter, Aq is the 
antenna amplitude, and Vg is the collision fre- 
quency of species s. Time is normalized to the lin- 
ear frequency, determined by the coUisionless gy - 
rokinetic dispersion relation ( Howes et al.ll2006l ). 
of a kinetic Alfven wave at the simulation do- 
main scale. Wo = 3.6kzoVA = 3.6a;^o- There- 
fore, the corresponding domain scale turn-around 



time is tq 



2tt/(jJ{) 



1.75a; 



AO- 



Collision fre- 



quencies are chosen to prevent build-up of small- 
scale structure in velocity space but remain small 
enough not to alter the weakly collisional dynam- 
ics: Vs <§; Wo is satisfied, so the simulation is 
weakly collisional. The expansion parameter e sets 
the parallel domain size, determined by assuming 
a critically balanced cascade driven isotropically 
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Fig. 1. — (Color online) Time averaged one- 
dimensional trace magnetic energy spectrum 
(solid) and the spectrum of the square of the cur- 
rent density (dashed). 



at the energy injection scale, kiPi = 10 ^, cho- 



sen to model solar wind turbulence (jHowes et al 
l2008ah . The antenna amplitude is chosen to sat- 
isfy critical balance at the domain scale, so the 
simulation represents critically balanced, strong 
turbulence. Prior analy sis of a similar simulation 
( TenBarge et al.ll2013al ) demonstrates a magnetic 
energy spectrum in excellent agreement with the 
empirical form found from a large statistical sam- 
ple of dissipation r ange measurements in the solar 
wind ( Alexandrova et al. 2012) . 



3. Magnetic Energy Spectrum and Cur- 
rent Density Spectrum 

First we demonstrate that the magnetic energy 
spectrum from our simulation is consistent with 
measurements in the solar wind, and then present 
the perpendicular wavenumber spectrum of the 
square of the current density, p[k±). In Figure [TJ 
we plot (solid) the average one-dimensional trace 
magnetic energy spectrum EBik±), where the to- 
tal magnetic energy Eg ' — J dk±EB{k±). The 
average is performed over the steady-state evo- 
lution of the system, 1.5ro < t < 4.1to; error 
bars represent the variance over the same inter- 
val. This spectrum is quantitatively consistent 
with a large sample of measurements of the dis- 
sipation ran ge magnetic energy spec trum in the 
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solar wind (jAlexandrova et al.M2012l ). suggesting 
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Fig. 2. — (Color online) Parallel current density, 
jz, for a perpendicular plane at t ~ 2A3tq, with 
different band-pass filters applied: (a) unfiltered, 
(b) 5 < k^pi < 21, (c) 21 < k±p^ < 84, and 
(d) k±pi > 84. Contours of the parallel vector 
potential A^ are shown in (a). 



that this simulation contains the essential phys- 
ical ingredients underlying turbulence in the so- 
lar wind. Since j = |j| = |(c/47r)V x B|, we ex- 
pect to find j ex kB, and therefore the scaling of 
P{k±) = j k±d(f)dk\\p{k) should satisfy the re- 
lation j'^{k±) oc kj^EB(k±), as confirmed by the 
plotted spectrum (dashed) of j"(fcj^). 

4. Current Sheet Formation 

In a weakly collisional plasma, magnetic recon- 
nection in the presence of a guide magnetic field 
is expected to develop an electron diffusion re - 
gion at the scale k±pe ~ 1 (JBirn fc Priestj|2006h . 
where a strong current sheet is expected to form. 
In Figure [21 we plot (a) the parallel current den- 
sity jz{x, y) from a perpendicular cut through the 
simulation domain. To explore the contribution to 
the current from different scales, we present spa- 
tially band-pass filtered data: (b) 5 < k^pi < 21, 
(c) 21 < k±pi < 84, and (d) k±p^ > 84. The 
large-scale currents visible in panels (a) and (b) 
are dominated by the driving, which generates up- 
ward and downward propagating kinetic Alfven 
waves with k±pi = 5. Panel (c), whose filter is ap- 
proximately centered on the electron gyroradius, 
0.49 < k±pe < 1.95, shows that this 3D gyroki- 
netic simulation indeed produces current sheets 



at the pe scale, consistent with such development 
in a wide range of plasma turbulence simulations. 
These electron-scale diffusion regions are highly 
intermittent, both spatially and temporally, with 
a typical lifetime r < O.Itq. The lack of significant 
current density in panel (d) shows that current 
sheets do not form at scales k±pe > 2, confirm- 
ing that our simulation has sufficient perpendicu- 
lar spatial resolution to capture fully the current 
sheet dynamics. 

Note that we are using the general definition of 
current sheet as a t angential disc o ntinuity in the 



magne tic field, e.g., iParken ()1994I ): iBirn fc Priest 
(J2006I ). This is also consistent with the defini- 



tion of current sheet used in recent solar wind 



literat u re. e.g..lVasauez et al. (|2007[ ) ; iGreco et al 



(|2009l ): IOsmanetal.1 (|2011al ). We do not 



at- 
tempt differentiate between current sheets asso- 
ciated with reconnecting magnetic flux and those 
arising from interfering Alfven waves. 

It is also important to note that the simulation 
is driven by injecting domain-scale waves, which 
generate strong turbulence. The current sheets 
form self-consistently from the cascade of wave- 
driven turbulence and are not seeded or other- 
wise initialized to form. Therefore, we conclude 
that electron-scale current sheets form as a natu- 
ral consequence of Alfvenic turbulence in this 3D 
gyrokinetic simulation. 

5. Current Sheets and Electron Heating 

We next examine the relative contribution 
to the measured electron heating from wave- 
particle interactions and dissipation in current 
shee ts. The analytical equations for plasma heat- 
ing ( Howes et al.ll2006[l have been implemented as 
a diagnostic in AstroGK. Boltzmann's H Theorem 
states that the entropy increase neces sary for irre- 



versib le heating requires collisions (jHowes et al. 



20061 ) ■ so the collisional heating (plus a small 
amount of numerical dissipation) is used to mea- 
sure the heating rate of each plasma species. Over 
the range of scales simulated, k±pi S [5, 105], lit- 
tle ion heating occurs, so we focus on the electron 
heating. The collisional electron heating is given 

by 
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dz 






L^ 



2i, 



tOe 
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/9/i, 



ekj^ 



V dt 
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[ dt 



cohJ 



where /lek^ = he(kx,ky, z,v\\,v±,t) is the non- 
Boltzmann portion of the perturbed electron dis- 
tribution function and -Foe i s the equilibriura 
electron distribution , function (JHowes et al.ll2006t 
2010r) . The collisional electron 



Numata et al. 



heating as a function of perpendicular wavenum- 
ber Qeik±) is computed by summing over annular 
rings in the perpendicular plane such that the to- 
tal heating is given by Qe — J dk±Qe{k±). Note 
that, in the weakly collisional limit, the heat- 
i ng rate is indepen dent of the collision frequency 
(JHowes et al.l 120061 ) . For the steady-state evolu- 
tion of the simulation over 1.5to < t < 4.1to, the 
heatin g diagnostics recover to tal power balance to 
< 2% (|TenBarge eraLll2013a[) . 

To estimate the contribution of current sheet 
dissipation to the heating rate, we compare the 
fraction of volume occupied by current sheets to 
the electron heating rate as a function of time. 
The volume filling fraction of current sheets is 
computed as the percentage of the volume with 
current density j > jth, with a chosen threshold 
jth = jmax/'i, where jmax is the maximum cur- 
rent density over all time and space in the sim- 
ulation. Note that varying this threshold alters 
the magnitude of the filling fraction but not the 
form of its variation with time. In Figure |3] a), we 
plot the volume filling fraction in percent (black 
dashed) and the electron collisional heating rate, 
discussed above, Qe (red solid) as a function of 
time — all quantities in the figure have been in- 
tegrated over the entire simulation domain. In 
Figure [3] b) is plotted the boxcar averaged (over 
At = 0.13to) power injected into the plasma by 
the Langevin antenna (magenta dotted) and the 
total energy of the turbulent fluctuations in the 
simulation including the magnetic field and kinetic 
energies Ekaw = Eb^^ + -Esy + Eke (blue dash- 
dotted), where Eke ~ X^s '^s'^OsMs/2 and u^ is 
the fluid velocity of each species. 

The cross correlations between the electron 
heating rate and filling fraction (black solid), an- 
tenna power (magenta dotted), and Ekaw (blue 
dash-dotted) are plotted in Figure [3] c). The in- 
jected antenna power and the electron heating 
rate are poorly correlated, (max(Corr((5, jant))) = 
0.52 ± 0.02 — the mean and variance are calcu- 




10" 



Fig. 3. — (Color online) a) Volume filling fraction 
of current sheets satisfying j > jmax/S in percent 
(black dashed) and the electron coUisional heating 
rate Qe (red solid), b) The injected antenna power 
(magenta dotted) and the total turbulent energy, 
Ex AW: (blue dash-dotted), c) Cross correlations 
between the electron heating rate and filling frac- 
tion (black solid), antenna power (magenta dot- 
ted), and Ekaw (blue dash-dotted). 



lated from the present simulation and five other 
identically prepared simulations employing dif- 
ferent random number seeds for the turbulent 
driving. Similarly, the electron heating rate and 
total turbulent energy are not well correlated, 
(max(Corr(Q,£;^AVK))) = 0.78 ± 0.04, suggest- 
ing that the heating rate is not simply a function 
of the driving or magnitude of turbulent energy. 
On the other hand, the electron heating rate is 
well correlated with the current sheet filling frac- 
tion, (max(Corr(Q,n/,n))) = 0.91 ± 0.04. The 
strong correlation suggests that dissipation asso- 
ciated with current sheets plays an important role 
in heating the electrons. 

6. Electron Heating by Scale 

In FigurelH we present a plot of the electron col- 
lisional heating rate by perpendicular wavenum- 
ber, Qe{k±) (solid), averaged over I.Stq < i < 
4.1 To, where error bars represent the variance over 
the interval. Note that the instantaneous shape 
of this heating curve is similar to the average, 
and does not change significantly as current sheets 
form and dissipate. This plot shows that the elec- 
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Fig. 4. — (Color online) Measured heating of 
the electrons from the simulation by perpendic- 
ular wavenumber, Qe{k±) (solid), an estimate of 
the electron heating based on linear wave-particle 
damping Qwp{k±) (dotted), and the Ohmic heat- 
ing rate Qn{k±) (dashed). 



tron heating is nearly constant over all scales, with 
about half of the total heating occurring at scales 
k±Pe < 1- The turn-down at k±pi > 105 is an 
artifact due to the diminishing number of Fourier 
modes in the corner beyond the fully resolved sim- 
ulation domain. 

As a function of k±, we may predict the colli- 
sionless damping of the turbulent fluctuations by 
resonant wave-particle interactions in our simula- 
tion using Qwp{k±) = 2"f{k±)EKAwik±), where 7 
is the linear kinetic damping rate of kinetic Alfven 
waves (dominated by electron Landau damping). 
This prediction for the wave-particle interaction 
heating rate requires integration over ku , where 
parallel is with respect to the local magnetic field 
and is typically determined via structure functions 
or wavelets. To avoid the complications of deter- 
mining the local magnetic field direction, we use 
frequency as a proxy for the parallel wave vec- 
tor since uj oc ku for k inetic Alfven waves— an ap - 
proach also taken in iTenBarge fc Howed ( 20121 ). 
This prediction for wave-particle damping, plot- 
ted in Figure m (dotted), admits no free param- 
eters, yet it agrees quite well with the measured 
coUisional electron heating (solid): the integrated, 
total predicted electron heating is within 4% of 
the coUisional heating diagnosed in AstroGK. The 



slight excess of wave-particle damping at 5 < 
k±pi < 40 and excess of electron collisional heat- 
ing at k±pi > 40 is consistent with the action of 



the e lectron entropy cascade (jSchekochihin et al. 
20091 ) ■ Through the entropy cascade, energy re- 



moved by electron Landau damping at 5 < k±pi < 
40 is expected to appear as collisional heating at 
higher wavenumbers, consistent with the simula- 
tion results in Figure ID 

We also compute the Ohmic heating rate Q,, = 
•qj^, whe re 77 = 0.38(47r)i/p,d^/c^ is the Spitzer re- 



sistivity (jSpitzer fc HarmI 119531 ). Vei = v^ is the 
electron-ion collision frequency, df. — c/ujpe is 
the electron inertial length, and Up,, is the elec- 
tron plasma frequency. The Ohmic heating rate 
is plotted (dashed) in Figure SI falling well be- 
low the measured heating rate and the predicted 



collisionless damping rate, Qr, ^ Qu 



Qe 



Theory predicts negligible Ohmic heating for a 
weakly collisional plasma, since electron-ion colli- 
sions are insufficient to significantly heat the elec- 
trons. Therefore, Ohmic dissipation of the current 
cannot account for the observed electron heating 
in the simulation, despite the strong correlation 
between current sheet filling fraction and heating 
rate. 

7. Discussion 

A puzzling aspect of these results is that, de- 
spite the clear correlation between the electron 
heating rate and the volume filling fraction of cur- 
rent sheets, the electron heating as a function of 
wavenumber is well predicted assuming that Lan- 
dau damping is entirely responsible for the elec- 
tron heating. This unexpected agreement raises 
the interesting possibility that the dissipation of 
the current sheets in the simulation occurs entirely 
via Landau damping. 

The solution to this puzzle lies in the rela- 
tionship between the current and magnetic field, 
namely j ex kB. This relationship implies that re- 
gions of strong current correspond to regions with 
enhanced small-scale magnetic structure since the 
current is weighted toward small spatial scales. 
Since the linear kinetic damping rate increases 
with wavenumber, regions with enhanced small- 
scale magnetic structure will also correspond to re- 
gions of enhanced wave-particle damping. There- 
fore, regions of strong current may also correspond 




Fig. 5. — (Color online) Quantities as in Figure [3] 
with a high-pass filter applied to remove modes 
with k±pi < 10, except the antenna current. 

to regions of enhanced wave-particle damping. 

To test this hypothesis, we apply a high-pass 
filter to the data presented in Figure [3] retaining 
only modes with k±pi > 10 — no filter is applied 
to the antenna current. The result of the filter- 
ing is plotted in Figure [5l where all of the cross 
correlations exceed 0.9, with the exception of the 
antenna power, which remains poorly correlated. 

The picture suggested by this simulation is 
one in which current sheets are self-consistently 
formed by the interaction of kinetic Alfven wave- 
like fluctuations, where each of the fluctuations is 
dissipating at its Landau damping rate. There- 
fore, current sheet formation and dissipation is 
dominated by the evolution of the Alfvenic tur- 
bulence, and current sheets correspond to sites 
of enhanced dissipation and heating, as sug- 
gested by recent ana lyses of solar wind turbulence 
(jOsman et allbOllallH ). 



8. Conclusion 

For a wave-driven, 3D gyrokinetic simulation of 
the dissipation range, k±pi G [5, 105], we find that 
current sheets self-consistently develop, in quali- 
tative agreement with fiuid simulations. The elec- 
tron heating rate is well correlated with the vol- 
ume filling fraction of current sheets, suggesting 
that the dissipation of current sheets plays an im- 
portant role in the heating of electrons. However, 
the electron heating rate as a function of scale 



is well predicted by assuming that all dissipation 
is due to collisionless damping of the turbulent 
fluctuations via Landau resonance with the elec- 
trons. In the weakly coUisional plasma, Ohmic 
dissipation of current sheets is negligible. Due to 
the relationship between the current and magnetic 
fleld, significant current highlights regions with 
enhanced small-scale magnetic structure, which 
will be collisionlessly damped at a rate greater 
than surrounding plasma. This suggests that cur- 
rent sheets may correspond spatially to locations 
of enhanced dissipation and heating, regardless 
of whether that dissipation is due to collisionless 
wave-particle interactions or active magnetic re- 
connection. 

The authors thank Homa Karimabadi for help- 
ful suggestions. This work was supported by 
NASA grant NNX10AC91G and NSF CAREER 
grant AGS-1054061. This research used resources 
of the Oak Ridge Leadership Computing Facility, 
supported by the Office of Science of the U.S. DOE 
under Contract No. DE-AC05-00OR22725. This 
research was supported by an allocation of ad- 
vanced computing resources provided by the NSF, 
partly performed on Kraken at the National Insti- 
tute for Computational Sciences. 

REFERENCES 

Abel, I. G., Barnes, M., Cowley, S. C, Dorland, 
W., & Schekochihin, A. A. 2008, Phys. Plas- 
mas, 15, 122509 

Alexandrova, O., Lacombe, C, Mangeney, A., 
Grappin, R., & Maksimovic, M. 2012, ApJ, 760, 
121 

Barnes, M., Abel, I. G., Dorland, W., et al. 2009, 
Phys. Plasmas, 16, 072107 

Birn, J., & Priest, E. R., eds. 2006, Rcconncc- 
tion of Magnetic Fields (Cambridge University 
Press) 

Borovsky, J. E., & Denton, M. H. 2011, ApJ, 739, 
L61 

Dmitruk, P., Matthaeus, W. H., & Seenu, N. 2004, 
ApJ, 617, 667 

Egedal, J., Daughton, W., Drake, J. F., Katz, N., 
& Le, A. 2009, Phys. Plasmas, 16, 050701 



Egedal, J., Daughton, W., & Le, A. 2012, Nature 
Physics, 8, 321 

Egedal, J., Fox, W., Katz, N., et al. 2008, J. Geo- 
phys. Res., 113, 12207 

Egedal, J., Le, A., Zhu, Y., et al. 2010, Geo- 
phys. Res. Lett., 37, 10102 

Frieman, E. A., & Chen, L. 1982, Phys. Fluids, 
25, 502 

Greco, A., Matthaeus, W. H., Scrvidio, S., Chuy- 
chai. P., & Dmitruk, P. 2009, ApJ, 691, LI 11 

Howes, G. G. 2008, Phys. Plasmas, 15, 055904 

Howes, G. G., Cowley, S. C, Dorland, W., et al. 
2006, ApJ, 651, 590 

— . 2008a, J. Geophys. Res., 113, A05103 

Howes, G. G., Dorland, W., Cowley, S. C, et al. 
2008b, Phys. Rev. Lett., 100, 065004 

Howes, G. G., TenBarge, J. M., & Dorland, W. 
2011a, Phys. Plasmas, 18, 102305 

Howes, G. G., TenBarge, J. M., Dorland, W., et al. 
2011b, Phys. Rev. Lett., 107, 035004 

Karimabadi, H., Roytershteyn, V., Wan, M., et al. 
2013, Phys. Plasmas, 20, 012303 

Le, A., Egedal, J., Daughton, W., et al. 2010, Geo- 
phys. Res. Lett., 37, 3106 

Lourciro, N. F., Schekochihin, A. A., & Zocco, A. 
2013, ArXiv e-prints 

Markovskn, S. A., & Vasquez, B. J. 2011, ApJ, 
739, 22 

Matthaeus, W. H., & Velh, M. 2011, 
Space Sci. Rev., 160, 145 

Numata, R., Howes, G. G., Tatsuno, T., Barnes, 
M., & Dorland, W. 2010, J. Comp. Phys., 229, 
9347 

Osman, K. T., Matthaeus, W. H., Greco, A., & 
Servidio, S. 2011a, ApJ, 727, Lll 

Osman, K. T., Matthaeus, W. H., Wan, M., & 
Rappazzo, A. F. 2011b, ArXiv c-prints 



Parashar, T. N., Shay, M. A., Cassak, P. A., & 
Matthaeus, W. H. 2009, Phys. Plasmas, 16, 
032310 

Parker, E. N. 1994, Spontaneous current sheets 
in magnetic fields : with applications to stellar 
x-rays. Vol. 1 (Oxford University Press) 

Pritchett, P. L. 2006, Geophys. Res. Lett., 33, 
13104 

Salem, C. S., Howes, G. G., Sundkvist, D., et al. 
2012, ApJ, 745, L9 

Schekochihin, A. A., Cowley, S. C., Dorland, W., 
et al. 2009, ApJS, 182, 310 

Servidio, S., Dmitruk, P., Greco, A., et al. 2011, 
Nonlin. Proc. Geophys., 18, 675 

Spitzer, L., & Harm, R. 1953, Phys. Rep., 89, 977 

TenBarge, J. M., & Howes, G. G. 2012, Phys. Plas- 
mas, 19, 055901 

TenBarge, J. M., Howes, G. G., & Dorland, W. 
2013a, ApJ, submitted 

TenBarge, J. M., Howes, G. G., Dorland, W., & 
Hammett, G. W. 2013b, J. Comp. Phys., in 
preparation 

Vasquez, B. J., Abramenko, V. I., Haggerty, D. K., 
& Smith, C. W. 2007, J. Geophys. Res., 112, 
11102 



This 2-coluimi preprint was prepared with the AAS lAT^rjX 
maeros v5.2. 



